The extension of radiative viscosity to superfluid matteiQ 
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The radiative viscosity of superfluid npe matter is studied, and it is found that to the lowest or- 
der of 5fi/T the ratio of radiative viscosity to bulk viscosity is the same as that of the normal matter. 
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As one of the most important transport coefficients, 
bulk viscosities of simple npe matter, of hyperon matter 
and even of quark matter, both in normal and superfluid 
states, have been extensively studied (H-[l8|. for more ref- 
erences see [l9| . 

In fact, the mechanical energy of density perturbations 
is not only dissipated to heat via bulk viscosity, but also is 
radiated away via neutrinos, this was first pointed out by 
Finiz and Wolf in 1968[l|. However, the damping mech- 
anism through neutrinos is ignored for several decades, 
until recently, Sa'd and Schaffner-Bielich [2(| named this 
mechanism the radiative viscosity, and found it is 1.5 
times larger than the bulk viscosity to all Urea processes 
in the lowest order of Sfi/T, both in nuclear matter and 
quark matter. Yang et al. [2l[ studied non-linear ef- 
fect of radiative viscosity of npe matter in neutron stars 
for both direct Urea process and modified Urea process, 
and found that non-linear effect will decrease the ratio of 
radiative viscosity to bulk viscosity from 1.5 to 0.5 (for 
direct Urea process) and 0.375 (for modified Urea pro- 
cess); which means that for small oscillations of neutron 
star the large fraction of oscillation energy is emitted as 
neutrinos, but for large enough ones bulk viscous dissi- 
pation dominates. 

It's well known that below certain critical tempera- 
ture nucleons in neutron star matter will be in superfluid 
states. The bulk viscosity of superfluid nucleon matter 
have been studied d, and it turns out that the 

superfluidity may strongly reduce bulk viscosity. This 
paper aims to study the radiative viscosity of superfluid 
nucleon matter, in other words, we will give the relation- 
ship between radiative viscosity and bulk viscosity in the 
superfluid case. 

Both the bulk viscosity and the radiative viscosity of 
npe matter are related to direct Urea process (n — > p+e+ 
77 e , p + e — > n+u e ) and modified Urea process (n + N — > 
p + N + e + V ei p + N + e — > n + N + v e ) in different 
conditions. As shown by [22|, the direct Urea process 
is allowed by the momentum conservation when pp n < 
Pf + PFc and for pure npe matter where pf = PF e , it 
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corresponds to n p /n > 1/9. This happens if the density 
is several times larger than the standard nuclear matter 
density po — 2.8 x I0 14 gcm~ 3 . In the following, we focus 
on the direct Urea process. 

Let us first recall the formulae of radiative viscosity of 
simple non-superfluid npe matter 0, H3| • Considering a 
periodic perturbation to the baryon number density 



n b (t) = n w + Re(5n b e lut ). 



(1) 



It will lead to a deviation from /3-equilibrium character- 
ized by 



Su = fl p + U e - fl n — 8 flp + 8 (l e - Sp n , 



(2) 



where fi n , fip and p,e are the chemical potentials of the 
neutrons, protons and electrons. Sfi can be expressed in 
terms of the variations of two independent variables 5n b 
and 6X„, 



Sft = C— + B6X P 
n b 



(3) 



where X p — n p /n b is the proton fraction and the coeffi- 
cient functions C and B are given by 



C = 



9 Hp 
p <9n B 



B = n b 
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(4) 
(5) 



To the leading order, the net reaction rate and the incre- 
ments of neutrino emissivity due to off-equilibrium could 
be written as 



r„ - Ttj 



-a. 



5h 



where \2U\2 



£loss — 5q^ 2 



Ao = ^e(T )0) , 



c 1071 (t n\ 

5o = l5T e(T ' 0) ' 



(6) 
(7) 

(8) 
(9) 
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£ = 5/i/T, and e(T, 0) is the neutrino emissivity in equi- 
librium 



e(T,0) = 3.3 x 1(T 14 f^J ' T 6 MeV 5 . (10) 

For a periodic process, the expansion and contraction 
of the system will induce not only the dissipation of os- 
cillation energy to heat, but also the loss of oscillation 
energy through an increasing of the neutrino emissivity. 
Bulk viscous coefficient C and radiative viscous coefficient 
1Z can be defined for the description of these dissipation 
mechanisms, respectively 



" 1 dt(V-vf 



(^diss ) 



(floss) = - / dt(V-vf, 



(11) 



(12) 



where v is the hydrodynamic velocity associated with 
the density oscillations, and r = 2it /uj is the oscillation 
period. Using the continuity equation, one obtains 



K = 2 <MS 3 (0 

where the energy dissipation is 



(£diss) = / (r„ - T-^Sfidt, 
Jo 



(13) 



(14) 



(15) 



and the neutrino emissivity caused by the oscillation is 

(floss) = / floss dt. (16) 
JO 

Here, we only present the results, for detailed calcula- 
tions see [201 



AC 2 



uo 2 + (XB/n b y 



From Eqs. ©and©, one gets 

n = l 



(17) 



(18) 



(19) 



This relation only holds for small oscillations. If the per- 
turbation amplitude is large enough, the non-linear effect 

must be taken into account and this simple relation is no 

i — i 

longer correct [21] . 



Now let us consider the effect of nucleon superfluidity 
on radiative viscosity. Assuming 



6/j, ' 



-'loss 



Apparently, 



n 



(20) 



(21) 



(22) 



In the following, we will show how to calculate S/X in 
the lowest order of £. 

In non-beta equilibrium, the net reaction rate of direct 
Urea process with nucleon superfluidity is 24 1 



-in 
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n / dQ o 



6(P f --p i )\M fi \ 2 l[p Fi 

3=1 



(2tt) 8 



and the total neutrino emissivity is[25|, [26[ 

3 



(23) 



47T 

+ 00 



rpKi 



n / 



3 = 1 



dx v x i v \J{x v - £, vj) + J{x v + £, Vj)], 



(24) 



where 



J(x,Vj) 



dx 1 dx 2 dx 3 f(z 1 )f(z 2 )f(x 3 ) 



xS(zi + z 2 + x 3 - x), 



(25) 



and the subscript j = 1,2,3 corresponds to n,p,e re- 
spectively. Vj = is the gap amplitude and Zj 



T 

x v = Iff is the dimen- 



(where j = 1,2), x 3 = 

sionless energy of the neutrino, and f(x) = (1 + e x )~ 1 is 
the Fermi-Dirac functions of nucleons and electrons, p Fj 
is the Fermi momentum and m* is the effective particle 
mass. dClj is the solid angle element in the direction of 
the particle momentum, and |M/i| 2 is the squared reac- 
tion amplitude. 

Here, we want to stress that z\ and z 2 in the above 
two formula carry all the information about nucleon su- 
perfluidity. Whether for neutron superfluidity or proton 
superfluidity, near the Fermi surface we have 



/i = sign(ry) a/(5 2 + if, 



(26) 
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where t] = vf(j> — Pf), and p-p are the Fermi velocity 
and Fermi momentum, respectively; and 5 2 — A 2 F(i)), 
A is the energy gap and describes the dependence 

of the gap on the angle •& between the quantization axis 
and the particle momentum. For different types of nu- 
cleon superfluidity, the expression of F($) is completely 
different, which can be seen in [24-26J. 

In the case of £ <C 1, to the lowest order we have 



J(x u — €,Vj)-J(x u + €,Vj) = -2£dJ{x v ,Vj)/dx v , (27) 



J(x„-£, Vj)+J(x v +^,Vj)-2J(x l/ , Vj) = £ 2 <9 2 J(ir„, Vj)/dxl, 

(28) 

then the net reaction rate is 
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5(P f -Pi)\M fi \ 2 l[p F . 
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x (~4£) / dx v x v J{x v ,Vj), 



(29) 



and the neutrino emissivity due to the departure from 
/3-equilibrium is 



floss = e(£,Uj) - e(0,Vj) 

3 
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x6£ y dx v x v J(x V: Vj). 
o 



(30) 



Note that, the leading order of Ar is the first order of £, 
while for £i OBS it is the second order of £. 

Unlike normal matter, Eqs. ([29)) and ([30f haven't exact 
analytical solutions. However, one can find that ATSfi is 
in proportion to £i oss , thus we can easily obtain 



S/X 



(31) 



which is the same as So/Xo. This means that the nucleon 
superfluidity doesn't change the value of lZ/(. 

In summary, we have studied the radiative viscous co- 
efficient of superfluid npe matter, and find that for di- 
rect Urea process, the ratio of radiative viscosity to bulk 
viscosity in the lowest order of £ yielded by Sa'd and 
Schaffner-Bielich [2(| for normal nucleons could be ex- 
tended to the superfluid case. In fact, it can be seen 
from our calculations that this is correct for modified 
Urea process, too. Thus, to the lowest order of £, the 
relation factor of 1 between the radiative viscosity and 
the bulk viscosity of npe matter is a generic one. 

The following are the extensive discussions to our re- 
sult: 

First, Although we have S/X = So/X in the low- 
est order, the viscosities in the superfluid case are dif- 
ferent from these in the non-superfluid case. The bulk 
and radiative viscous coefficients could be easily ex- 
pressed as £ = Co-Rr and 1Z = TZoR e , where £o and 
IZq arc the viscosities in the non-superfluid case, i?r = 
Ar(£,^)/Ar(£) and R e = 4>ss(£, ^)M° SS (0 are the 
reduction factors caused by superfluidity. Of course, our 
result shows that to the lowest order of £, R e equals to 
i?r (note that Rr has been calculated numerically by 
Haensel et al. 0,0). 

Second, our result satisfies the general relationship 
d£\ oss / ddpb — 3(r,7— r„), which has been found by Flores- 
Tulian and Reisenegger [27| . They found that this rela- 
tionship holds both in the case of normal nucleons and 
in the case of superfluid nucleons, and both in linear case 
and for precise solutions. 

Finally, our result based on the fact that to the low- 
est order of £, the existence of superfluidity don't change 
the relationship between the neutrino emissivity due to 
the departure from /3-equilibrium and energy dissipation 
in the form of heat. Nevertheless, it is not correct for 
higher order calculations. As shown in [2l],[23j], for nor- 
mal matter Ar and £i OS s can be solved analytically as 
polynomials of £, while for superfluid matter they must 
be solved numerically. As a result, we expect that in 
non-linear regime, the ratio of radiative viscosity to bulk 
viscosity of superfluid matter no longer equals to that 
of normal matter. However, the calculation in the non- 
linear case is far more complicated and we will consider 
it in our further study. 
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